Three-band s± Eliashberg theory and the superconducting gaps of iron pnictides 
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The experimental critical temperatures and gap values of the superconducting pnictides of both 
the 1111 and 122 families can be simultaneously reproduced within the Eliashberg theory by using 
a three-band model where the dominant role is played by interband interactions and the order 
parameter undergoes a sign reversal between hole and electron bands (s±-wave symmetry). High 
values of the electron-boson coupling constants and small typical boson energies (in agreement with 
experiments) are necessary to obtain the values of all the gaps and to correctly reproduce their 
temperature dependence. 
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The recently discovered Fe-based pnictide supercon- 
ductors [TJ [2j [3] have aroused great interest in the scien- 
tific community. They have indeed shown that high T c 
superconductivity does not uniquely belong to cuprates 
but can take place in Cu-free systems as well. Never- 
theless, as in cuprates, superconductivity occurs upon 
charge doping of a magnetic parent compound above a 
certain critical value. However, important differences ex- 
ist: the parent compound in cuprates is a Mott insulator 
with localized charge carriers and a strong Coulomb re- 
pulsion between electrons; in the pnictides, on the other 
hand, it is a bad metal and shows a tetragonal to or- 
thorhombic structural transition below w 140 K, followed 
by an antiferromagnetic (AF) spin-density-wave (SDW) 
order [4 . Charge doping gives rise to superconductiv- 
ity and, at the same time, inhibits the occurrence of 
both the static magnetic order and the structural transi- 
tion. The Fermi surface consists of two or three hole-like 
sheets around Y and two electron-like sheets around M. 
Up to now, the most intensively studied systems are the 
1111 compounds, ReFeAsOi_ x F x (Re = La, Sm, Nd, Pr, 
etc.) and especially the 122 ones, hole- or electron-doped 
AFe2As2 (A = Ba, Sr, Ca). The huge amount of exper- 
imental work already done in 122 compounds is due to 
the availability of rather big high-quality single crystals. 

Most of the present research effort is spent clarifying 
the microscopic pairing mechanism responsible for super- 
conductivity. The conventional phonon-mediated cou- 
pling mechanism cannot explain the observed high T c 
within standard Migdal-Eliashberg theory and the inclu- 
sion of multiband effects increases T c only marginally [5] . 
On the other hand, the magnetic nature of the parent 
compound seems to favor a magnetic origin of supercon- 
ductivity and a coupling mechanism based on nesting- 
related AF spin fluctuations has been proposed [6]. It 
predicts an interband sign reversal of the order parameter 
between different sheets of the Fermi surface (s± symme- 
try). The number, amplitude and symmetry of the super- 
conducting energy gaps are indeed fundamental physical 
quantities that any microscopic model of superconductiv- 
ity has to account for. Experiments with powerful tech- 



niques such as ARPES, point-contact spectroscopy, STM 
etc., have been carried out to study the superconducting 
gaps in pnictides (for a review see [7J). Although results 
are sometimes in disagreement with each other, a multi- 
gap scenario is emerging with evidence for rather high 
gap ratios, A1/A2 ~ 2 — 3 [7 . A two-band BCS model 
cannot account either for the amplitude of the experi- 
mental gaps and for their ratio. Three-band BCS models 
have been investigated [8j |9j [10] which can reproduce the 
experimental gap ratio but not the exact experimental 
gap values. In this regard a reliable study has to be car- 
ried out within the framework of the Eliashberg theory 
for strong coupling superconductors [11] , due to the pos- 
sible high values of the coupling constants necessary to 
explain the experimental data. 

By using this strong-coupling approach, we show here 
that the superconducting iron pnictides represent a case 
of dominant negative interband-channel superconductiv- 
ity (s±-wave symmetry) with high values of the electron- 
boson coupling constants and small typical boson ener- 
gies. Furthermore we prove that a small contribution of 
intraband coupling does not affect significantly the ob- 
tained results. The model is compared with the results 
of two representative experiments in 122 [12 and 1111 
[T3] compounds proving to be able to reproduce fairly 
well the values and the whole temperature dependence 
of the superconducting energy gaps. 
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FIG. 1: Schematic drawing of the multiband model used in 
this work. The two hole bands (1 and 2) are centered around 
the T point, while the equivalent electron band (3) around 
the M point of the reduced Brillouin zone. 
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As a starting point we can model the electronic struc- 
ture of pnictides by using a three-band model (FigfTJ 
with two hole bands (1 and 2) and one equivalent elec- 
tron band (3) [8 j. The s-wave order parameters of the 
hole bands have opposite sign with respect to that of the 
electron one [6] . Intraband coupling could be provided by 
phonons while interband coupling by antiferromagnetic 
spin fluctuations. In a one-band system spin fluctuations 
(SF) are always pair breaking but in a multiband one 
the interband term can contribute to increase the critical 
temperature. Indeed, in the multiband Eliashberg equa- 
tions (EE) the SF term in the intraband channel has 
positive sign for the renormalization functions Zi and 
negative sign for the superconducting order parameters 
A^ thus leading to a strong reduction of T c . However, if 
we consider negative interband contributions in the A^ 
equations, the final result can be an increase of the crit- 
ical temperature [T4] . 

Let us consider the generalization of the Eliashberg 
theory [11] for multiband systems, that has already been 
used with success to study the MgB 2 superconductor [T5] . 
To obtain the gaps and the critical temperature within 
the s-wave, three-band Eliashberg equations one has to 
solve six coupled integral equations for the gaps Ai(iuj n ) 
and the renormalization functions Zi(iu) n ), where i is a 
band index that ranges between 1 and 3 (see Fig[T]) and 
uo n are the Matsubara frequencies. For completeness we 
included in the equations the non-magnetic and magnetic 
impurity scattering rates in the Born approximation, Yf- 
and T%: 

UJnZ^lUJn) = C0 n + ^(1^ + T% )N? \lU n ) + (1) 
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where Ag(zu; n , zcj m ) = A^(^ n ,^ m ) + A*?(iu n ,iu m ), 

A£j(iu n ,iu m ) = A?j(iu) n ,iu) m ) - A*?(iu n ,iu m ). 
6 is the Heaviside function and oj c is a cut- 
off energy. In particular, A P ^ ,sp (icj n ,icj m ) = 

J +o ° dna^FP h ' a P(Q)/[{u; n - com) 2 + ^ 2 ], where ph 
means "phonon" and sp "spin fluctuations". Fi- 
nally, N^(iu m ) = A j (uj m )/^u 2 n + A 2 (iuj m ) and 

Nf(lU m ) = UJ m j ^uj 2 m + A|(zcj m ). 

In principle the solution of the three-band EE shown in 
eqsjl] and |2] requires a huge number of input parameters: 

i) nine electron-phonon spectral functions, a 2 jF ph (fl); 

ii) nine electron-SF spectral functions, a 2 ,F sp (Q); hi) 



nine elements of the Coulomb pseudopotential matrix, 
H*j{uj c )\ iv) nine non-magnetic (T^-) and nine paramag- 
netic (r^ ) impurity scattering rates. 
It is obvious that a practical solution of these equations 
requires a drastic reduction in the number of free param- 
eters of the model. On the other hand, from the work 
of Mazin et al. [8 we know that: i) X^ » X P ^ « 
i.e. phonons mainly provide intraband coupling but the 
total electron-phonon coupling constant S^A^ 1 should be 
very small [5], ii) X^ p >> X^ p « 0, i.e. SF mainly pro- 
vide interband coupling. We include these features in the 
most simple three-band model by posing: X^ = X P ^ = 0, 
X^ p = and /i^(a; c ) = fi*j(u; c ) = 0. In addition, we set 
r£ = r# = 0ineqs. (1) and (2). 

Under these approximations, the electron-boson 
coupling-constant matrix is then [8]: 

X 31 V! \ 
A 32 ^2 
A31 A32 J 

where v x = JVi (0)/JV 3 (0) , v 2 = N 2 (0)/N 3 (0) and A^(0) 
is the normal density of states at the Fermi level for the 
i-band (i = 1, 2, 3 according to Fig.l). 

We initially solved the EE on the imaginary axis to 
calculate the critical temperature and, by means of the 
technique of the Pade approximants, to obtain the low- 
temperature value of the gaps. In presence of a strong 
coupling interaction or of impurities, however, the value 
of Ai(icu n =o) obtained by solving the imaginary- axis EE 
can be very different from the value of A^ obtained from 
the real-axis EE [16 . Therefore, in order to determine 
the exact temperature dependence of the gaps, we then 
solved the three-band EE in the real-axis formulation. 

We tried to reproduce the critical temperature and the 
gap values in two representative cases: i) the 122 com- 
pound Ba .6Ko.4Fe 2 As2 with T c = 37 K where ARPES 
measurements gave Ai(0) « 12.1 meV, A 2 (0) w 5.5 
meV and A 3 (0) « 12.8 meV [12]; ii) the 1111 com- 
pound SmFeAsOo.sFo.2 with T c = 52 K where from 
point-contact spectroscopy measurements we obtained 
Ai(0) = 18 ± 3 meV and A 2 (0) = 6.15 ± 0.45 meV Q3]. 

Inelastic neutron-scattering experiments suggest that 
the typical boson energy possibly responsible for super- 
conductivity ranges roughly between 10 and 30 meV 
[T7] . In our numerical simulations we used spectral func- 
tions with Lorentzian shape, i.e. a 2 -F(Q) = Cij[L(Q + 
ilij.Yi^-L^-Q.i^Yij)] where L(n±n ij ,Y ij ) = [(ft± 
^j) 2 + (^j) 2 ] -1 ' are the normalization constants nec- 
essary to obtain the proper values of A^-, while f^j and 
Yij are the peak energies and half- widths, respectively. 
In all our calculations we always set f^j = Qq, with Qq 
ranging between 5 and 35 meV and Y^ = 2 meV. The 
cut-off energy is uo c — 12 • ^0 an d the maximum quasi- 
particle energy is uo max = 16 • ^o- 

In the 122 case (T c = 37 K) we know that v\ = 1 and 
v 2 — 2 [8] while in the 1111 case (T c = 52 K) we have 
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FIG. 2: Full symbols, left axis: Calculated gap values at T—2 
K for Bao.6Ko.4Fe2As2 (upper panel) and SmFeAsOo.sFo.2 
(lower panel) as function of typical boson energy Qo- Open 
symbols, right axis: Electron-boson coupling constants, A31 
and A32 as function of Qo • The inset shows the spectral func- 
tion used in this model in the case 0»j = 10 meV. 

v\ = 0.4 and V2 = 0.5 [18 . Once the energy of the boson 
peak, ^0 is set, only two free parameters are left in the 
model: A31 and A32. 

By properly selecting the values of these parameters 
it is relatively easy to obtain the experimental values of 
the critical temperature and of the small gap. It is more 
difficult to reproduce the values of the large gaps of band 
1 and 3 since, due to the high 2Ai ) 3//c#X' c ratio (of the 
order of 8-9), high values of the coupling constants and 
small boson energies are required. Fig |2] shows the values 
of the calculated gaps (full symbols, left axis) as function 
of the boson peak energy, Qq. The corresponding values 
of A31 and A32, chosen in order to reproduce the values of 
T c and of the small gap, A2, are also shown in the figure 
(open symbols, right axis). In both materials, only when 
Qo < 10 meV the value of the large gaps correspond to 
the experimental data. Indeed, when Qq increases, the 
values of Ai and A3 strongly decrease. As a consequence, 
a rather small energy of the boson peak together with a 
very strong coupling (particularly in the 3-1 channel) is 
needed in order to obtain the experimental T c and the 
correct gap values. In this regard, it is worth noticing 
that the absolute values of the large gaps cannot be re- 
produced in a interband-only, two-band Eliashberg model 
[T9] as well as within a three-band BCS model. In the 
latter case it is only possible to obtain a ratio of the gaps 
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FIG. 3: Calculated temperature dependence of the gaps 
for Bao.6Ko.4Fe2As2 (T c = 37 K, upper panel) and for 
SmFeAsOo.sFo.2 (T c = 52 K, lower panel): Ai(T) (red solid 
line), A 2 (T) (blue dashed line) and A 3 (T) (green dash-dot 
line). Symbols are experimental data from ref. [12] (upper 
panel) and ref. [13] (lower panel). 

close to the experimental one [9] [10] . 

Another important result of the model is the temper- 
ature dependence of the gaps. Figure [3] shows this de- 
pendence for the experimental gaps (symbols) together 
with the theoretical A^(T) curves obtained by the three- 
band Eliashberg model (lines) for Bao.6Ko.4Fe2As2 (up- 
per panel) and SmFeAsOo.sFo.2 (lower panel). The pa- 
rameters used for the 122 compound are ^0 = 10 meV, 
A31 = 4.267 and A32 = 0.569; for the 1111 compound we 
used Q = 10 mev, X 31 = 14.520 and A 32 = 1.708. The 
experimental temperature dependence of the gaps shown 
in the upper panel is rather unusual with the gaps slightly 
decreasing with increasing temperature until they sud- 
denly drop close to T c . The theory reproduces very well 
this behavior, which is possible only in a very strong cou- 
pling regime [16] . The different temperature dependence 
observed in the lower panel of Fig. [3] results from a com- 
plex non-linear dependence of A$ vs. T curves on A31. 
Further details will be given in a forthcoming paper. 



A31 = 3.866 


A 32 = 0.471 


fiij 0, \u 0.4 


Ai = 10.30 meV 


A 2 = 5.62 meV 


|A 3 | = 10.24 meV 



TABLE I: The effect of a small contribution of intraband cou- 
pling, Xu = 0.4 for the case of Bao.6Ko.4Fe2As2 at T — 2 K. 
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A 3 i = 2.730 


A 32 = 0.758 


rfj = 0.1, Xu = 0.4 


Ai = 7.49 meV 


A 2 = 5.72 meV 


|A 3 | = 7.98 meV 
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TABLE II: The effect of the Coulomb interaction, fi*j for the 
case shown in Table [I] 

We also tested the effect into the model of a small 
intraband coupling (possibly of phonon origin). In the 
case of Bao.6Ko.4Fe2 As2 we used A^ = 0.4 since we know 
indeed that this coupling cannot be very high [5]. It 
might be thought that this term can sensibly contribute 
to increase the gap values but, as can be seen in Table 
[IJ this is not the case: the gap values only show a slight 
increase (of the order of 1%). 

The effect of Coulomb interaction was also investigated 
for the case shown in Table [J where a weak intraband 
coupling is included. We chose fi*j = 0.1 and, as ex- 
pected, we found that the intraband Coulomb pseudopo- 
tential has a negligible effect while the interband one [T4] 
strongly contributes to raise T c and reduces in a consid- 
erable way the value of A31: in this case, as shown in 
Table [TTJ it is only possible to obtain the correct value 
of the small gap. As a consequence, this result seems to 
exclude a strong interband Coulomb interaction in these 
compounds. 

Finally we have also examined, for SmFeAsO0.sF0.25 
the case of a spectral function with two peaks at energies 
Qi and ^2- In the upper panel of Fig. [4] the two boson 
energies are f^i = 10 meV and ^2 — 20 meV, while in the 
lower panel we have Qi = 10 meV and ^2 = 30 meV. The 
gap values (left axis) and the coupling constants, A31 and 
A32 (right axis) are plotted as a function of the weight 
of the low-energy peak, w p . As somehow expected, when 
the weight of this peak is larger, (w p = 0.75) the gaps Ai 
and A 3 are larger and close to the experimental ones but 
the coupling constants, A31 and A32 strongly increase. 

In conclusion, we have shown that the newly discov- 
ered iron pnictides very likely represent a case of dom- 
inant negative interband-channel pairing superconduc- 
tivity where an electron-boson coupling, such as the 
electron-SF one, can became a fundamental ingredient 
to increase T c in a multiband strong-coupling picture. In 
particular, the present results prove that a simple three- 
band model in strong-coupling regime can reproduce in 
a quantitative way the experimental T c and the energy 
gaps of the pnictide superconductors with only two free 
parameters, A31 and A32, provided that the typical ener- 
gies of the spectral functions are of the order of 10 meV 
and the coupling constants are very high (A31 > 4). 

We thank I.I. Mazin and E. Cappelluti for useful dis- 
cussions. 
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FIG. 4: The calculated gaps and electron-boson coupling con- 
stants A31 and A32 in the 1111 case as function of the weight 
w p of the low-energy peak Q±. Q 2 = 20 meV (upper panel) 
and 30 meV (lower panel). 
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